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ABSTRACT 



In order to obtain the dynamical structure of the tropical 
atmosphere in equatorial latitudes a short range prediction experiment 
is proposed. The initial field contains the intertropical convergence 
zone and associated disturbances over western Pacific Ocean during 
March 1965* A complete initial state is constructed using a consistent 
balance system of equations . The procedure involves construction of 
pressure, temperature and vertical motion distributions starting from 
an observed rotational part of the wind field. It is shown that this 
procedure, whose validity assumes a small Rossby number, does not yield 
a realistic structure of the vertical motion. A short range prediction 
with a multi-level prediction model yields some interesting solutions 
in the vicinity of the intertropical convergence zone. During the 
first 18 hours of prediction an adjustment of the motion and the mass 
field ensues with gravity inertia oscillations. A detailed discussion 
of some dynamical aspects at 2b hours is presented. An important 
feature of the model is a parameterization of cumulus convection as a 
function of large-scale moisture convergence. The role of cumulus 
scale heating in the vicinity of the intertropical convergence zone 
is explored by carrying out experiments with and without diabatic 



heating. 



List of symbols 



a 


fraction (or percent) of synoptic scale area covered by 
convective clouds 


C D 


drag coefficient 


C 

P 


specific heat of air at constant pressure 


e 


vapor pressure of air. Subscript to refers to water surface 


E v 


evaporation from sea surface 


f ’ f o 


Coriolis parameter, and a mean value of the Coriolis parameter 
( respectively) 


F , F 
x y 


friction force per unit mass of air 


F s 


flux of sensible heat from ocean to the atmosphere 


CD 

& 


Diffusion of heat and moisture, respectively 


g 


acceleration of gravity 


H, M 


heat and moisture sources and sinks, respectively 


V h l 


sensible and latent heat (respectively) per unit mass of air 


h 


terrain height above sea level 



, j|, Ik unit vectors along zonal, meridional and vertical directions 



i 


net moisture convergence in a unit vertical column extending 
from 900 millibars to the top of the atmosphere (100 mb) 


ITCZ 


intertropical convergence zone 




eddy kinetic energy 


K z 


kinetic energy of mean zonal flow 


l 


number of grid points along x-axis 


L 


latent heat of vaporization of air 


m 


map projection factor 


P, P Q 


pressure and a standard pressure (= 1000 mb), respectively 


P E 


eddy potential energy 



eddy potential energy 



Q. 

Q, 



R 



q 

Ro 

t, e 

T 

a 



T S’ 



U, V, 0) 

V, \V 



0 ) 



C 



0) 

T) 

mn 

x,y,p>t 



potential energy of mean zonal flow 
any arbitrary variable 

amount of moisture required to replace an ambient moisture 
sounding by a local moist adiabat 

amount of moisture required for condensation warming to 

replace an ambient temperature sounding by a local moist adiabat 

gas constant 

specific humidity 

Rossby number 

temperature and potential temperature of air, respectively 
temperature of air at anemometer level, assumed same as ^qqq 
in this study 

temperature and specific humidity (respectively) of a parcel 
raised vertically above the 900 millibar surface with no 
lateral mixing 
temperature of water 

zonal, meridional and vertical (p) velocity 

total horizontal wind speed and wind vector, respectively 

average vertical velocity over region occupied by convective 

clouds inside a synoptic-scale elementary area 

average vertical velocity over a region not occupied by 

convective clouds inside a synoptic-scale elementary area 

synoptic scale vertical velocity 

coefficient of a 9-point local surface 

independent variables 



Z 



height of constant pressure surface 



0 

e 


equivalent potential temperature 


U , < 


horizontal and vertical eddy diffusion coefficients 


T , T 

X y 


surface frictional stresses along the x and y axes 


Ap 


vertical grid distance 




geopotential height of pressure surface 




zonal mean geopotential on a boundary wall 


$ 


latitude 


(i, J) 


a reference grid point on the mercator map 


P 


density of air 


* *1000 


the subscript 1000 refers to the value at the 1000 millibar 
surface 


<l», c 


streamfunction and relative vorticity 


% 


non-divergent part of the wind 


0 


static stability parameter, = - — 




absolute vorticity 


6 


the beta parameter, - 


n 


RT 

Exner pressure, = — 


X 


divergent part of the wind field, a velocity potential 


( " ) 


bar refers to a zonal average 


( )' 


primes at a symbol refer to departures from a zonal average 



1. Introduction, 



In the middle latitudes the formation of an extratropical storm 
is generally regarded as a baroclinic instability problem. The primary 
energy source that leads to large kinetic energy of extratropical vor- 
tices is the available potential energy. Phillips (1956) discussed the 
energetics of this class of disturbances in his numerical general circu- 
lation experiment. 

Lack of understanding of the types of instabilities in the tropics 
has, so far, prevented the use of well defined numerical experiments to 
simulate the development of disturbances that are frequently portrayed 
from observations on weather maps. 

Some success in the understanding of the formation and maintenance 
of a tropical storm (or hurricane) has been achieved through the works of 
Ooyama (1963), Charney and Eliassen (196*0, Kuo ( 1965 ) and several others. 
The primary reason for the success of this work has been due to a recogni- 
tion of what is called a conditional instability of the unsaturated 
tropical air mass. Unlike the baroclinic instability problem where dis- 
turbance with horizontal wavelengths of the order of 1000 or more kilo- 
meters have a maximum growth rate, for the tropical conditional instabil- 
ity the maximum growth rate is found for horizontal wavelengths of 5 to 
10 kilometers in scale (Kuo, 1961). This latter scale is now recognized 
as the cumulus scale. Most active tropical disturbances contain this 
cumulus scale in the form of organized bands of convection. There is a 
certain cooperation of the large and small scale the nature of which is 
a central problem in tropical meteorology. The entire belt of tropics is 
frequently conditionally unstable, namely 
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(0e is the equivalent potential temperature and p is pressure.) This is 
observed from an examination of atmospheric soundings below the 500 mb 
surface. However, convective scale as manifested by satellite and air- 
craft observations of clouds, is present only over portions of the tropics. 
Thus, there is evidently a supression of the convective motion over large 
areas of the tropics. This might be controlled by the large scale atmos- 
pheric dynamics. 

A mere examination of the streamline features on a weather map does 
not, in general, give indications of where one should expect convective 
motion to be present or absent. Several tropical meteorologists present 
large quantities of weather maps showing merely the field of wind direction. 
The atmosphere being considerably more complicated, not many questions can 
be thus either asked or answered. 

Because of large data voids in the tropics the age of discovery of 
new phenomena is not yet over. As an example of this one may cite the 
discovery of the climatological mean troughs of the upper tropical tropo- 
sphere over the Atlantic and Pacific Oceans. These were deduced from com- 
mercial aircraft reports of upper winds during recent observational 
studies by Dean (1956) and Sadler (1967). These papers do not deal with 
either the three-dimensional structure or the dynamics of these climato- 
logical features. There is thus no wide recognition of the importance of 
these studies. The data is sufficient to portray cyclonic or anticyclonic 
streamline features on a constant pressure surface, but it is generally 
insufficient for study of the dynamical structure. Several climatological 
features of this kind must exist over the remaining vast uncharted oceanic 
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regions of the tropics. Through use of satellite cloud photographs it 
might be possible to locate regions where such disturbances are present. 

This paper presents some features of the tropical circulations 
for one day, March 1, 19^5, as part of an investigation which is being 
continued over a longer period of time. The study deals with the dy- 
namics of the intertropical convergence zone. The approach is to con- 
struct a consistent three-dimensional distribution of motion, mass, 
temperature and moisture fields through use of diagnostic and prognostic 
numerical models. The large-scale horizontal wind speed in the tropics 
is of the order of 10 meters/sec while vertical motion, w , is of the 
order l cm/sec (Baumhefner, 1968). The latter may be expressed by a 
relation of the form 

_ . 

a) = a a) + (l-a) a) (l) 

c 

where , a l is fraction of the synoptic scale area covered by convective 
clouds, and u) are respectively vertical motions in the cloud and out- 
side the cloud in a unit area. Theories of parameterization of the small 
scale convective motion as a function of a larger scale motions enable 
us to determine the parameter f a T subject to a number of approximations, 
to can be estimated from solution of balance equations and adequation or 
through numerical integration of primitive equations. 

No description of the three-dimensional structure of a weather 
disturbance can be complete without a distribution of the motion field 
in the three dimensions. Any investigation of energetics, development 
and instability relies heavily on the distribution of vertical motion. 

For phenomenon of small Rossby number (Ho < 1), a measure of the 
vertical motion may be obtained from a solution of the quasi-geostrophic 
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to-equation or a more general balance a)-equation (Krishnamurti , 1968a). 

In investigations of phenomenon of low latitudes (like the intertropical 
convergence zone) where Ro *,1, the validity of the quasi -geostrophic or 
balance m-equation is questionable. Primitive equation forecast models 
contain no approximation (except the quasi-static approximation, if the 
vertical equation of motion is a hydrostatic balance) and hence a short 
range forecast (using a quasi-geostrophic or balanced initial state) 
might yield a reasonable solution of the three-dimensional motion field 
which would be valid for large Rossby number. Generally, such an unbal- 
anced initial state is expected to produce gravity inertia oscillations. 
These are damped oscillations and there is an ensuing adjustment of the 
dependent variables into a large Rossby number disturbance. Some map 
features of the initial state are lost as both the wind and the pressure 
field undergo a degree of adjustment. The process of adjustment of wind 
and pressure in low latitudes in the presence of the smaller convective 
scale is a very complicated and important theoretical problem. Simple 
formulations of this problem by Rosenthal (1965) and Matsuno (1966a) are 
important contributions. 

We shall show in this paper, that the numerical solution of the 
primitive equations (starting from solutions of a diagnostic balance 
model) undergoes an adjustment, the theoretical nature of which is not 
quite clear, but the process is certainly of interest. It should be 
noted that the adjustment and the damping of the gravity inertial oscilla- 
tion is brought about by the diffusion of momentum and potential tempera- 
ture. The diffusion coefficients are selected somewhat arbitrarily in 
these studies*, the choice was made by trial and error. In our study after 
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a period of approximately 12 hours these oscillations become very small 
in amplitude. We are thus able to present a description of the three- 
dimensional distribution of mass, moisture and motion field for large 
Rossby number phenomena through short range numerical prediction of 
atmospheric variables. Two questions that come up in the light of these 
comments regarding the oscillations are: 

i) How good is the initial data for a description of the large 
scale distribution of the dependent variables in regions near the equa- 
tor or in the tropics? 

ii) How much of the initial map features is lost during this 
process of adjustment? 

We are concerned here with a study of the intertropical converg- 
ence zone. There is a considerable amount of data in the western Pacific 
Ocean to describe the large-scale flow both at the surface level and up 
to the 200 millibar surface. This is borne out by a number of descript- 
ive studies in this region (Yanai, 196l). This data is supplemented by 
a number of aircraft reports (AIREPS) that enable a better analysis of 
the 700 and 250 millibar charts. The initial analysis at the lowest 
levels contains features like a convergence line (the ITCZ). At the upper 
level features like warm high pressure cells are found to the north and 
south of the low level convergence line. The anticyclonic outflow from 
these upper level high pressure cells (size lOOOnds of kilometer scale) 
is an easterly wind maximum over the equatorial latitudes and lies 
directly above the surface level convergence asymptotes. The damped 
oscillation of the prediction does not destroy these large-scale features. 
Hence, a description of the predicted fields may be considered a description 
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of the large-scale dynamical structure of the ITC. This is the principal 
motivation for carrying out this investigation. 

2. A primitive Equation Model 

(a) Advective scheme and basic equations: 

The main features of the advective scheme are essentially the 
same as those proposed by the author, Krishnamurti , 1962. Similar quasi- 
Lagrangian schemes have been proposed by Leith ( 196 U ) and Okland ( 19^5 ) - 
Leith carried out long-term numerical integrations in a general circula- 
tions model while Okiand proposed its use in short-term integrations for 
a two-level model. No non-linear advective terms (\V • VQ) appear explic- 
itly in this formulation and thus no obvious non-linear computational 
instabilities are present. In a recent review of various problems gener- 
ally encountered in numerical weather prediction, Elsasser (1968) points 
out that (jkland’s Lagrangian advective scheme was one of the important 
factors in the successful prediction of cyclogenesis over short time 
periods of the order of two days. 

The main outline of advective scheme will be discussed here. 
Consider the basic system of the primitive equations expressed in the 
form 
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This may be regarded as a closed system for the seven variables u, v, 0 , 

3 z 

^1000’ “ 3p ° r alternativel y a closed system for u, v, 0, q, z and 
a), with appropriate boundary conditions. is a quasi-Lagrangian opera- 



tor given by 



_D_ 

Dt 




m 




+ 




( 9 ) 



Here m is a map factor, in our case for a mercator projection expressed 
by the relation 

m = sin 85 cosec (90 - $) (10) 

where $ is the latitude, m = 1 at $ = 15 °N and 15 °S 5 at the equator 
m = 0 . 996 . Equations (2) and ( 3 ) are the momentum equations, 

(U) is a form of the thermodynamic energy equation. Equation ( 5 ) is a 
statement for conservations of water vapor. Equation (6) is a tendency 
equation of the pressure frame where the role of orography, h, is in- 
cluded. (7) and (8) are, respectively, the continuity and hydrostatic 
equations. Heat sources and sinks, H, and the moisture sources and 
sinks, M, are discussed later in the paper. F^ and F , F Q , and F^ repre- 
sent frictional and diffusive transfers of momentum, heat and moisture, 



and are also discussed later. 
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The principle of the quasi-Lagrangian advective scheme is to 
locate the past position (P) of a parcel on a pressure surface during 
each time step, the parcel arrives at a grid point 0 (I, J) at the end 
of each time step, see figure 1. During the course of movement from 
P to Q non-conservative changes take place which are incorporated in 
the advected properties. 

The iterative procedure utilizes a local surface over 9 points 
of figure 1 defined by the relation 




2 2 



l l 



n X® Y n 
mn 



( 11 ) 



The nine coefficients ri^ are predetermined for each point and do not 
vary during the course of the integration. X and Y are distances in 
the zonal and meridional directions between the central point I, J and 
the other 8 points of figure 1. 

Distance (PO.) = tyx 2 + Y 2 . 

Point P is determined by a predictor corrector successive ap- 
proximation procedure, the first guess is defined by the relations 

X * - ' 






( 12 ) 



Y ' " At 

The variables u, v, w, z, 0, q. A, B, C, D and E are determined at the 
point P using relation (ll) and the values of X and Y defined in 
equation (12). 

Dote that At is a time step: the subscript t-At is the previous time 
period* and t would be the current time for which a prediction is 
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desired. The 


first guess prediction of the 


dependent variables is 


carried out according to the relations, 
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(7) and (8). 










The second guess defines a location 


of the point P through the 


relations 
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, Z, 0, and q are interpolated 


at point P (time t-At) 


and the prediction at point Q is carried out using the relations 
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-< P >t-At * B(9) t 4t 
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Z(Q) t and u)(Q) t are again determined from equations (7) and ( 8 ). 

This two-step predictor corrector is similar to Matsuno’s (1966b) 
time differencing scheme for Eulerian and spectral equations ( see also 
Lilly (1965)). Terms containing vertical advection, pressure gradients, 
Coriolis force, friction, heating and other non-conservative features on 
the right hand side of equations (2) through (7) are incorporated in A, 

B, C, D and E. A time step At = 10 minutes is used in the present study. 

b) Domain of integration, distribution of variables: 

The domain of integration extends from 10U°E longitude to 128 °W 
longitude over the Pacific. Northern and southern limits of the domain 
extend from 28°N to 28°S. The region between 152°W to 128°W is used to 
prescribe a cyclic continuity in the east-west direction. This is 
carried out by a quadratic interpolation of the initial analyzed data 
such that , 

Q(10U°E) = Q(128°U) (16) 

where Q is any analyzed variable. It must be noted that the initial 
state over the entire domain is a consistent balanced hydrostatic set, 
cyclic continuity merely removes east-west boundaries in both the dia- 
gnostic and prognostic studies presented here. 

In a recent comparative study of several finite difference 
analogs for primitive equations, Grammeltvedt ( 1968 ) states that if in- 
tegrations are performed for more than three days then more than 15 grid 
points per wavelength are needed to describe with accuracy the movement 
and development of the shortest wave which initially is carrying a sig- 
nificant part of the energy. This is an important result for numerical 
weather prediction. This would imply that for studies of development and 
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movement of a disturbance in the tropics, whose wavelength is of the 
order of 1000 km, a mesh size of the order of 6 0 km would be required. 

A fine mesh model will be widely used, when faster computers with larger 
core storage become available. Phase speed errors, through reduction of 
truncation errors, can be reduced appreciably by this procedure. In the 
present study we use a mesh size of 2° latitude by 2° longitude and its 
limitations are thus obvious. In certain large-scale general circulation 
models, mesh sizes of the order of 5° latitude are being used. 



horizontal; however, there is a staggering of dependent variables in the 
vertical. It is a five-level model and staggering permits all nine levels 
to be used. This vertical structure is primarily dictated by the nature 
of the dynamical and thermodynamic equations, see figure 2. 

c) Boundary conditions for the primitive equation model 

As mentioned earlier, there are no boundaries in the zonal direc- 
tion in this study. Our basic grid extends from 28°N to 28°S. There is, 
however, an implicit north and south pole in the model. At the north and 
south poles we assume at each time step 



In the present study there is no staggering of variables in the 



z iooo (North pole) = l Z 1000 (28 ° N)/£ 

X 



6 (North pole) = £ 6 (28°N)/£, 



X 



q (North pole) = £ q (28°N)/£ 



x 




(17) 




x 



0 (South pole) = J 0 (28°S)/£ 



x 



q (South pole) = J q (28 °S)/£ 



x 
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u, v and to at the two poles are assumed zero during the integra- 
tion. A linear interpolation between pole and 28° is used to calculate 
the values at 30°N and 30°S during each time step. These values were 
used as boundary conditions during the course of the integration. This 
scheme was found to be stable and values of the vertical motions, the 
most sensitive parameter of this model, do not become unreasonable. This 
boundary condition does not in any way represent the behavior of the real 
atmosphere at 30°M and 30°S. 



d) The development in the primitive equation model 

In quasi-geostrophic models development is generally discussed in 
terms of processes producing individual changes of vorticity via the vor- 
ticity equation. These changes are primarily attributable to the pres- 
ence of fields of divergence and convergence. In primitive equation 
models surface development may be described via the follbwing equations. 

rlOOO 

< ia > 

100 



“lOOO 



_D_ 

Dt 



( Z 1000 - h ) 



“1000 pg 6 1000 



(18b) 



If h = 0 (no terrain) and if the vertically integrated divergence 
is zero, then 



and 



“lOOO 0 



DZ 1000 

Dt 



0 



(19) 



( 20 ) 



i.e., parcels moving on the 1000-millibar surface will conserve their 
heights and no deepening or filling of surface systems is possible. 
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Up and downslope motions over large-scale terrain features permit de- 
velopment of surface systems such as upslope anticyclogenesis and a 

lee cyclogenesis. Net convergence leads to a sinking motion ( a) ^000 > ^ 

DZ 

and — > O) > while the converse is true of net divergence. One 
Dt 

would think that since the vertically integrated divergence is a small 
difference between large numbers (at each level) the results of such 
vertical integration will be very sensitive and only spurious results 
can be obtained. This, however, definitely is not the case. A bal- 
anced initial state starts out with minor gravity inertia oscillations 
which damp rapidly with time and very small and reasonable magnitudes 
of ^2000 are produced by the time integrations, see figure 12. 

This has been verified by nearly all of the research workers who have 
integrated the primitive equations. A synoptic meteorologist finds 
it very difficult to understand this feature of the primitive equation 
models when he knows that he has only a limited success in the use of 
the so-called tendency equation. The apparent reason for this lies in 
the fact that a primitive equation model generates a dynamically con- 
sistent momentum field with sufficient accuracy such that the net di- 
vergence and tendencies are small, as is the case for the real atmosphere. 
Direct observations do not yield small net divergences or tendencies 
because winds usually are in error. A detailed, dynamic explanation of 
this process of adjustment would be of considerable synoptic interest. 

In the vicinity of mountains w^oOO acquires magnitudes very close 
to that given by the lower boundary conditions of the m-equation , i.e. 

^1000 = ~ ^ p 1000 W 1000 ‘ This, again, is somewhat surprising because 

in the primitive equation model, terrain h modifies ^000 v ^ a a ra ^her 
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complicated cycle of operations. 

h modifies Z^qqq via equation (6) 

^1000 ^ e ^ erm: *- nes Z via equation (8) 

Z determines u, v via equations (2) and (3) 
u, v determines ^000 v ^ a equation (7) 

The reason for this must again be that a dynamically consistent 
primitive balance is indeed obtained when the initial state is given by 
the balance equations. The nature of such a primitive balance in the 
vicinity of the equator, where the Rossby number is greater than unity, 
has not been explored in theory as yet. 

3. Initialization 

A primitive equation model is currently being used in operational 
weather prediction at the National Meteorological Center in the United 
States (Shuman and Hovermale, 1968). The initial wind field is obtained 
from a solution of the so-called linear balance equation. This relation- 
ship between the wind field and the pressure field does not permit any 
divergence or vertical motions at the initial time. A similar initial 
state was proposed for the Bushby and Timpson's ( 1967 ) model by Benwell 
and Bretherton (1967) in a recent paper. It is pointed out by the lat- 
ter authors that this procedure suppresses inertio-gravitational oscilla- 
tions. These ideas are by no means new, Charney (1955), Hinklemann (1951) 
and Smagorinsky et al (1965) have proposed the use of non-divergent ini- 
tial states for integration of primitive equations. Phillips (i960), on 
the other hand, showed that an initial divergence equal to that prescribed 
by quasi-geostrophic dynamics would be required to suppress the unwanted 
gravity inertio-oscillations . 
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In this paper we propose a rather detailed diagnostic balance 
model to define a consistent initial state. Thus, initially both diver- 
gence and vertical motions are present. A part of our interest will be 
to describe the three-dimensional synoptic structure and the dynamics of 
the motion, mass, thermal and moisture variables at the initial state 
from balance equations. This study departs somewhat from earlier study 
on initialization by the author (Krishnamurti , 1968a); here we propose 
an initialization of the dependent variables across the equator in a 
tropical belt. This poses seme problems if the streamfunction were to 
be obtained from observed geopotential heights, but it is a relatively 
simple matter if the streamfunction is deduced from the observed wind 
field. The initialization procedure is as follows: 

(a) Streamfunction ^ 

The nondivergent part of the wind is obtained from a solution 
of the equations 

v 2 ip = c (21) 

and \V = |k x Vip (22) 



where £ is the relative vorticity of the observed wind. Boundary condi- 
tions for i|> are obtained by the procedure outlined by Hawkins and Rosen- 
thal (1965) (version 2 of their work). This particular choice of boundary 
condition permits no net mass flux from the domain; however, i|> is permitted 
to vary along the boundaries . 

(b) Initial distribution of geopotential heights 

The distribution of geopotential ,<J> = gZ, is obtained by solving 
the balance equation 



VI = 



, . m + 2J (ft , ft) 



(23) 



16 



Since i() (x, y, p) is known the right hand side is defined in the domain 
and the solution of this elliptic equation is a straightforward problem. 
At the boundaries we assume <j> to be given by the relation 

<{> B = 4> + f Q (\p - (23a) 

4> b is the value of <j> at the northern or southern walls at 28°N or 28°S. 

$ is the zonal mean value of <p obtained from observations of geopotential 

height of the northern or southern wall. 

-k -1 

fg = 10 sec at northern wall, 

-if -i 

and fg = - 10 sec at southern wall. 

\p is the value of the streamfunction at any point along the wall and 
ip the zonal mean value along that wall. 

The initial distribution of potential temperature is given by 




Potential temperatures computed by this procedure were found to depart 
by as much as 3 degrees in the middle latitudes from observations. 
Maximum departures were of the order of 5 degrees at lower levels. 

This difficulty could have been overcome by using observed tempera- 
tures directly in the model; this, however, would cause the sacrifice 
of the initial balance relationship of the variables. It is not 
a priori clear how serious an effect this would have had on the initial 
amplitudes of gravity inertio-oscillations . 

( c ) Divergent part of the wind, X, and vertical velocity 

A detailed analysis of the system of equations required to solve 
for X and oj is presented elsewhere by the author (Krishnamurti , 1968a). 
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Here we shall merely present the equations and their boundary conditions. 
Omega equation 



_2 ^ .2 3 2 oj 

V oto + f — r- 

3p 



f '4 J( ' M a ) + "’ 2 J< *- 6) + k 



3r 3 t 
3x 3 y 



- 2 — J fli lih f JL frn 2 y) 

2 3t 3p J ^3x’ 3y-* f *- C7 X ' 



3p 



r-p ’ 2 " * f ,7 (" £ ^ 



* r F» < 5 “ • v IS - f £ f' x • 



ap 



3p 



- n v 2 ( 7 x.ve}- 6 Ai|l 



(25) 



Continuity equation 



v 2 x = 

3p 



(26) 



Vorticity equation 



V 2 If = - J(*.s a ) * w 



VC + c V X 
a a 



3 „2, „ „ 3if> 3 

- o> ^ V * - Vm • V - - g - 



3F 3F 

y x 

3x 3y 



(27) 



Boundary conditions 



At p = 100 mb, a) = 0 



At southern wall, y = y^, to = X = |^- = 0 



At northern wall , y = y , uj = X = = 0 

C. O t 



At p = 1000 mb , to = - g 



1000 



R 0 



1000 



j(i|), h) - VX • Vh 



■] 
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The solutions for w, X and are obtained by a successive approxi- 

<3t 

mation procedure. This involves three-dimensional relaxation for m and 
two-dimensional relaxation for X and . These three equations are 
coupled. Conditions for ellipticity of these equations are discussed by 
Krishnamurti (1968a). 

and are frictional forces per unit mass of air, along the 
x and y axis:, H is net diabatic heating per unit mass of air. These are 
discussed in some detail in a separate section of this paper. The same 
forms of friction and heating are retained in the initialization and pre- 
diction equations. 

The total initial wind is given by 



\V = 




(28) 



Initial distribution of moisture is a scalar hand analysis of 
relative humidity at standard meteorological levels, this is converted 
into fields of specific humidity q at the 900, TOO, 500 and 300 millibars 
surfaces. The closed system of initial dependent variables are u, v, m, 
Z, 6 and q, figure 2 shows the vertical staggering of variables. 



I 4 . Diabatic heating and parameterization of cumulus scale 

In the present model diabatic heating per unit mass of air is 
split into two components. 

H = H s + H l (29) 

where H is the transfer of sensible heat from the ocean: 
is the release of latent heat. 

Ue have not considered radiative effects. 
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(a) Sensible heat flux : 

A simple flux formula (Jacobs, 1951), was used to define this 
quantity. 



Fg = 0.0329 V (T w - T fl ) (30) 

The unit of flux, F , is mb meter sec”\ V being the wind speed 
in meters/sec at anemometer level and (t - T ) the temperature differ- 
ence between water at the ocean's surface and air at the anemometer 
level. The constant of proportionality has dimensions mb/°C. The rate 
of heating Hg may be expressed by the relation 



H 



S 




(31) 



Potential temperature 6 is a dependent variable and it appears 
at the 900, 700, 500 and 300 millibar surfaces. It is assumed that 
sensible heat transfer modifies the potential temperature of the air at 
the 900 millibar surface, and the flux, Fg, is assumed zero at 800 mb 
and above. 

We express sensible heat at 900 millibar surface by the formula 

H S = fefT °- 0329 V < T „ - T a> 

= 0.00323 V (T - T ) (32) 

w a 

2 -3 

Here H has units meter sec . If the air-sea temperature difference 

is of the order of 1°C and the wind speed is of the order 10 mps, then 

a change of one degree in potential temperature may he expected over a 

2U-hour period. This would require a representative analysis of the 

initial distribution of (T - T ) over the oceans. Because of general 

w a 
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lack of data only a rather smooth representation of this field is pos- 
sible at this time. Radiative cooling is a serious emission of the 
present study because it is also of the order of a degree centigrade 
per day. 

(b) Latent heat 

H t may again be divided into two forms of latent heat, stable 
and conditionally unstable. The former may be important in regions of 
dynamic ascent of moist and absolutely stable air. The latter may be 
important if there are regions where certain prescribed conditions for 
the evolution of cumulus scale convective elements are met. For the 
present we have omitted the stable form of latent heat as the tropical 
air mass is generally known to be conditionally unstable in the lower 
troposphere and our interest in this study is an equatorial, tropical 
belt . 



The procedure for parameterization of the cumulus scale as a 
function of the large-scale motion with the mutual modification of the 
two scales is carried out along the lines of Kuo's (1965) formulation 
for a tropical storm. 

Heat is released if the following criteria are met: 

i) _ — < 0 at any level above 1000 millibars. This measures 

3 P 

the moist conditional instability. 

ii) Net moisture convergence, I > 0, in vertical columns. The 
column extends from the top of the friction layer to the top of atmos- 
phere (100 mb). 



The heating function, H , is defined by the relation: 

L 



H 



L 




- T) 



( 33 ) 
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Where (T - T) is the temperature difference between cloud element and 

D 



tropical environment r, is the specific heat of air at constant pres- 
sure; At is a cloud time scale, where in a time integration problem we 
use At = the time step of the primitive equation model (Kuo, 1965); and 
"a" is a measure of the fraction of the grid scale area covered by act- 
ive convective elements. 

The calculation of "a" involves determining the following three 
integrals . 

Net moisture convergence above friction layer. 

W 900 q 900 



I = * 

G 



r P T 



q V dp - 



900 
■ P T 



g 



On 



gk (a -s - «> * 



(3U) 



(35) 



900 



* C ^ • 



(36) 



and 



a = I/(Q 1 + Q 2 ) 



(37) 



T and q are respectively the temperature and specific humidity 

O D 

of a parcel raised dry adiabatically from 900 millibar surface to the 
lifting condensation level and then moist adiabatically to a level p^ 
where intersects the environmental temperature. A detailed analysis 
of this parameterization procedure and aspects of calculation are dis- 
cussed elsewhere by the author (Krishnamurti , 1968b). In this earlier 
study diagnostic calculations of vertical motions and various friction 
layer parameters were evaluated. It was shown that convective banded 
motions as revealed by airborne radar and that calculation of vertically 
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integrated moisture convergence and the field of the parameter "a" 
exhibit reasonable correlations. In the initialization and prediction 
phase of the present study the heating function is evaluated according 
to the flow chart, figure 3, which is self explanatory. 

Large scale dynamical prediction of tropical disturbances includ- 
ing a formal parameterization of the cumulus convection have been carried 
out by Manabe and Smagor insky ( 1967 ). It would be pertinent to review 
their parameterization procedure. In their general circulation model 
moist convective heating and condensation takes place when the lapse 
rate tends to become super-moist adiabatic and relative humidity tends 
to exceed 100 per cent. During this situation a convective adjustment 
is forced upon the numerical solution in such a way as to make the re- 
lative humidity and the equivalent potential temperature uniform in the 
unstable layer. Furthermore this adjustment is carried out in such a 
way as to keep the sum of potential, internal, and latent energy invari- 
ant. One of the obvious limitations of this procedure is that no con- 
vective process on small scale is permitted if relative humidity < 100 
per cent. Synoptic experience in the tropics is quite to the contrary 
(e.g. Riehl, 19^5; Baumhefner, 1968 ). Manabe and Smagorinsky ( 1967 ) 
discuss a rather interesting tropical disturbance that forms in their 
long-term general circulation experiment. Its large-scale structure is 
quite similar to that of a weak tropical storm, although a large inner 
core of this disturbance is very moist. The parameterization procedure 
that has been adopted here permits convective heat and moisture flux in 
the vertical direction even when the relative humidity is less than 100 
per cent. However, this scheme presently suffers frem the lack of a 
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built in check on convective adjustment which would control vertical 
soundings and force them to look like their atmospheric counterparts. 

Such a built-in convective adjustment is considered essential for long- 
term general circulations experiments. In this experiment calculations 
were extended by time steps to U8 hours and no severe problems such as 
supersaturations or super-moist adiabatic instabilities were found. 

Kuo’s scheme, used in this study, behaves roughly in the following manner 

A conditionally unstable region initially experiencing a net 
convergence of moisture flux, I > 0, is modified in a matter of a few 
hours so that large-scale values of T approach Tg, and q approaches qg. 
This causes "a" to become smaller and the region to become sufficiently 
moist that further mass flux into the region begins to bring in relative- 
ly drier air and I changes sign to I < 0. When this happens we cease to 
parameterize the convective scale. 

(c) Evaporation and precipitation 

Conservation of water vapor is given by equation (5). Where M 
represents the net sources and sinks of water per unit mass of air. 

Rohwer (l93l) expresses evaporation from the ocean surface, E , 

v 

by the relation: 

E v = (2.6 ♦ 7.7 V 10Q0 ) (0.98 e u - e 1000 >/(2**- * 3600 >. (38) 

2 

The units of are gm/(cm sec); V^ooo I s 1000 millibar wind 
in meters/sec; e^ is the saturation vapor pressure in millibars at the 
temperature of the water’s surface; and e is the vapor pressure of 
air at the 1000 millibar surface. No attempt was made to vary the rough- 
ness constant for variable wind speeds over the oceans. Over land areas 
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was set equal to zero. 



Webb (i960) has critically reviewed the use of such an empirical 
formula for estimating evaporation and states that studies over Lake 
Eucumbe show that error would be within 10 per cent. 

Precipitation was assumed to be entirely convective and was evalu- 
ated at each layer using Kuo f s parameterization procedure, figure 3. 



K 



C a (T - T) 
P S 

LAt 



( 39 ) 



K refers to the vertical level and P has the units sec - "*". 

K 

In our calculations the atmosphere is divided into four layers, 
1000-800, 800 - 600 , 600-U00 and kOO to 200 millibars. Evaporation is only 
permitted in the lowest layer, and M is accordingly defined by the 
relation : 



M 



E C a (T_ - T) 

V p o 

g Ap ” LAt 



( 1 * 0 ) 



with M = - 



C a (T_ - T) 

_E § 



LAt 

of the lowest layer. 



in the other three layers. Ap is the thickness 



5 . Surface friction and diffusion terms of model 

These terms, F^ and F^, in the momentum equations contain two 
parts in this model: 

i) Subgrid scale lateral and vertical diffusion of momentum and 
ii) effect of surface wind stress. 

This formulation is almost similar to that used by Mintz ( 196 U ) in his 
general circulation studies. 



We write : 
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_ 2 ^2 3^u 3t x 

F x = pm V u + < — + g — 

3p 



and 



2 2 3^v ^ T y 

F y = 7 v + < “a + 8 jf 

3p 



(Ul) 



( 1 * 2 ) 



U is a lateral eddy diffusion coefficient whose value in the present 

U 2-1 

experiment is 5 x 10 meter sec and k is a vertical eddy diffusion 

-2 2 -I 

coefficient whose value in the present experiment is 2 x 10 mb sec 

The primary role of sub-grid scale lateral and vertical diffusion of 

momentum is to remove energy from the smallest permitted scales. It 

cannot represent physical subgrid scale processes in the horizontal or 

vertical, such as mesoscale horizontal jets or cumulus convection. It 

is found that absence of vertical diffusion k gives rise to rather large 

horizontal wind speeds at the lowest levels within a ^8-hour forecast 

period. Calculations did not appear to be critically dependent on the 

choice of the lateral coefficient. 

Surface frictional stresses t and t are defined in terms of 

x y 

a surface roughness constant applied at the lower level, p = 1000 mb. 



2 

= C n pi/u' + v 



u 



( 1 * 3 ) 



T 

y 






u 2 + 



2 

v v 



(kk) 



Cp is the drag coefficient, whose value in the present experiment is 
2.5 X 10 -3 . 

Surface stresses prescribed in this manner provide a crude 
boundary layer for our model between the 1000 and 900 millibar surfaces. 
Vertical motions at the 900 millibars are influenced significantly by 
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the presence of this effect. Frictional convergence of mass enhances 
low-level convergence of moisture and the latter prescribes a percent 
area of convective cloud distribution in our model. A discussion of 
these aspects is presented with the results of the numerical experiment. 

6. A description of the Initial State: March 1, 1965 > 1200 Z. 

It is not possible to describe all of the scalar fields that 
define the initial state. A selection of some of the pertinent fields 
is presented here. 

a) Streamlines of the observed wind field at 1000 and 200 millibar 
surfaces . Figures 4a and 4b, respectively, portray the streamlines 
(solid lines) and isotachs (dashed lines) at these two levels. 

The following synoptic features of the intertropical conver- 
gence zone may be noted at the lower level. An asymptote of conver- 
gence extends from 152°W to about l44°E near 5° north latitude. The 
speed field along this asymptote is quite variable. Wind speed ranges 
from 5 to 25 knots over most of the map. There appears to be no simple 
structure for the isotach field in the vicinity of this asymptote of 
convergence. Wind speed maxima can be found along and north of this 
asymptote. 

There is a belt of easterlies over most of the northern hemi- 
spheric portion of the analysis. An asymptote of convergence near 
108°E is shown as a part of a cyclonic eddy at 3°N near l40°E. In the 
southern hemispheric tropics interesting features are the two cyclonic 
disturbances at 1T°S, 110°E and at 15°S, 129°E. The speed field is 
shown to be somewhat organized in these disturbances. Over the rest of 
the Southern Hemisphere there are no well marked centers. 



27 



At 200 mb there are two regions of westerly wind maxima, one 
in the northwestern portion of the map over eastern Asia, and another 
over central Australia in the Southern Hemisphere. Over the equatorial 
latitudes a broad easterly current with speeds near 1+0 knots prevails . 

Four anticyclonic circulation centers (two in each hemisphere) consti- 
tute an important part of the upper-level flow. There are several 
other transient centers that are portrayed in the streamline analysis 
of this chart. 

b) Moisture analysis at 900 m illib ar surface : Moisture analysis was 

carried out by first using all available stations at the surface level. 

A satellite nephanalysis (valid at approximately this time) was used as 
an underlay and the surface streamline features were used to complete 
the subjective distribution of relative humidity at 1000, 850, 700 and 
500 millibar surfaces. The upper levels (850, 700 and 500) were arrived 
at by subjective examination of cloud cover, upper level data and the 
surface analysis. The relative humidity does not exceed saturation 
values on the large scale at any point. The analysis is very subjective 
and the available information is merely a guide. The distribution of 
the moisture variable is very poorly defined, hence it is somewhat unre- 
presentative. The initial values of specific humidity at 900 (interpo- 
lated between 1000, 850 millibars) are shown in figure 5(a). A moisture 
maximum is found north of Australia in the equatorial Pacific Ocean. 

It has values of specific humidity > l8 gm/kgm. Another moisture maxima 
center with values close to 20 gm/kgm lies northeast of Australia and 
just south of the equator. Over the continent of Australia air is 
fairly dry, values being close to 5 gm/kgm. Over most of the equatorial 
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tropics values are near 15 gm/kgm. There is a gradual decrease of 
moisture between l6°N and 28°N to around 6 gm/kgm near 28°N. 

Because of the uncertainty in the analysis of the moisture vari- 
able a prediction experiment was also carried out with uniform initial 
relative humidity at each pressure surface. This experiment is dis- 
cussed separately below. 

c) Temperature distribution in the middle troposphere : 500 millibar 

temperatures are shown in figure 5(b). Thermal structure over the 
central Pacific is characterized by two nearly isothermal warm centers 
with temperature around 269°K. These warm regions define the warm core 
structure of the anticyclones that extends to the 200 millibar surface; 
see figure Mb). 

In the northern and southern portions of the maps two well- 
marked baroclinic zones are present . They define regions where wind 
maxima are found at the 200 millibar surface, figure 4(b). The active 
portion of the intertropical convergence zone lies in the vicinity of 
3°N and lU0°E. In the middle troposphere mesoscale temperature anomalies 
are likely to be present over this region. The smooth large-scale iso- 
thermal patterns could be interpreted either as a lack of smaller 
scale data or an efficient lateral and vertical mixing process. This 
could be possible if vertical momentum exchange in the convective re- 
gion reduces the vertical wind shear (Gray, 1967). 

Sea surface temperatures were obtained from the monthly mean 
charts produced by U.S. Bureau of Commercial Fisheries (1962) and 
Hoyskoi Atlas (1953) and were kept fixed during the experiment. Ter- 
rain height distributions were obtained from the tabulations of Berkofsky 
and Bertoni (i960). 
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d) Satellite cloud cover : During this period (March 1-3, 1965) TIROS X 

made a few passes and relevant pictures of cloud cover. Figure 6 (a and 
b) shows the cloud cover from approximately 130°E to l80°E and between 
20°N and 20°S. The intertropical convergence zone is very well marked 
with bright clouds near North between l60°E and l8o°E on March 1. 

It is somewhat more extended and brighter on the 2nd and again well 
marked at about the same latitude on the 3rd. In the vicinity of 170°E 
and 15°S a well-marked cloud band oriented northwest-southeast is observ- 
able on all three days. On the 1st and 2nd a bright band of cloud cover 
is also persistent in the region 130°E to ll*0°E and 10°S to 20°S. This 
is over north central and north western Australia. The three above-men- 
tioned regions of intense cloud cover are found in the vicinity of 
cyclonic surface systems and asymptotes of low-level convergence, figure 
U(a). Over other regions the amount of cloudiness is variable during 
the three-day period. The picture time in each case is within 6 hours 
of the synoptic time of analysis. 

No attempt will be made to relate each cloud feature observed on 
these TIROS passes with results of computations. Such an attempt would 
be futile at the present state of our work. Some reference will be made 
to the observed cloud cover, at least in its larger-scale context, in 
the following sections. 

Me shall next present some of the results of the initialization 
from the balance equations. 

e) Initial streamfunction and vorticity fields : Results are presented 

for the 200 and the 1000 millibars surfaces, figures 7(a and b), respec- 
tively. The solid isopleths in these figures show the streamfunction 
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and the dashed isopleths show the distribution of absolute vorticity 
(units 10 ^ sec ^). 

The streamfunction distribution at 200 mb contains one large 
anticyclonic circulation center in the northern Pacific Ocean and a 
center over northern Australia. Strong easterlies are found over the 
equatorial latitudes between these centers. At higher latitudes an upper 

troughs in the westerlies contain values of absolute vorticity ^1 x 10 

-1 -4 -1 

sec in the Northern Hemisphere and ^ - 1 x 10 sec in the Southern 

Hemisphere. The zero isopleth of absolute vorticity meanders across the 
equator. At the 1000 millibar surface a maximum value of absolute vor- 
ticity is found southeast of the Truk Islands in the northern Pacific 
Ocean. This is associated with a cyclonic disturbance along the western 
end of the asymptote of convergence. Observed wind directions and the 
streamfunction are nearly parallel at the 200 millibar surface, while 
they are not nearly so in the vicinity of the asymptote of convergence 
at the 1000 millibar surface. A stronger divergent component of wind 
was present at the lower level in the vicinity of the ITCZ. The cyclonic 
disturbances over the subtropical latitudes of the Southern Hemisphere 
are indicated by the large negative values of absolute vorticity. 
f ) Vertical Motions on top of the friction layer and calculated cloud 
cover : Figure 8 shows the isopleths of the initial balanced vertical 

motions x 10 ^ mb/sec at the 900 millibar surface (solid lines), while 
the calculated percent area of convective cloud cover is shown by dashed 
lines. The shaded region has some fraction of the synoptic scale grid 
(2° latitude x 2° longitude) covered by convective clouds, and the un- 
shaded region is free of convective clouds. The subtropical disturb- 
ances of the Southern Hemisphere exhibit marked convective cloudiness. 
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Maximum active cloud amounts range from 0.5 to 1 per cent in these dis- 
turbances. Balanced vertical motion in these disturbances has a maximum 
value os- 60 x 10~^ mb/sec which is o^l cm/sec. The principal region 
of interest in this study is the surface asymptote of convergence. As 
stated earlier, the validity of the balance vertical motions near the 
equator are questionable because Rossby number ^1.0. 

A wide region of maximum rising motion (- 60 x 10 ^ mb/sec) with 
maximum cloud amounts ^1 per cent is found in the vicinity of the ITCZ. 
Satellite cloud cover shows a narrow band of clouds, figure 6* (a and b), 
along U°N near 170°E and a region with widespread convective activity 
between 1^0°E and 150°E. A region of sinking motion extends from 136°E 
to 152°E along the equator, this implies that air reaching the equator 
from the Southern Hemisphere is sinking at the rate of 1 cra/sec before 
it begins to ascend in the vicinity of the asymptote of convergence in 
the Northern Hemisphere. These calculations of cloud cover and the in- 
tensity of vertical motion appear somewhat unrealistic in the vicinity 
of the equator. Short-range prediction with this initial distribution 
of mass, motion and moisture variables reveals some interesting changes 
in the vicinity of the equator; this is discussed in the following 
section. 

7. A discussion of the time evolution of various fields 

Some of the results of the calculations will be discussed in 
this section. 

a) 1000 and 200 millibar wind field at 2k hour s : 

Figures 9 (a and b), respectively, show these fields. The plotted 
winds are the observations at 2k hours, and the isopleths are the predicted 
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are extremely interesting. The wind field at 1000 millibars exhibits 
an east-west oriented asymptote of convergence. Convective activity 
is rather well marked along this asymptote of convergence. We felt that 
the forecasts after 30 hours were not very reliable, mainly because of 
boundary errors. The major asymptote of convergence and associated con- 
vective cloud band exhibits an equatorward movement. Convective activ- 
ity weakens somewhat as the band moves from U°N at 2k hours to near the 
equator. At 30 hours this system strengthens again and is located at 
5°S at U8 hours of forecast time. Meanwhile a new asymptote of converg- 
ence forms at 12 to l6°N latitudes. Our calculations show that rather 
large changes can occur in the tropics during a ^8-hour period. This is 
not surprising because satellite photographs, figure 6 (a and b) , indi- 
cate marked changes in cloud cover. Manabe and Smagorinsky (1967) also 
mention that tropical disturbances have a rather variable structure 
during 1 to 3 day periods. 

c ) Vertical motions at 900 millibars (10 ^ millibars/sec) 

During the course of the time integration other relevant fields 
like convergence of mass between 1000 and 900 millibars, convergence of 
moisture between 1000 and 900 millibars and the net convergence of 
moisture between 1000 and 100 millibars were printed out at 6-hour in- 
tervals. It is true, that numerical models cannot be expected to reveal 
any more physics than what is already stated by the basic equations of 
our models. We did not impose any very direct relationships in our 
basic equations between the cloud cover computations and the various 
parameterization features listed above. We observed the following re- 
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Asymptotes of convergence, percent cloud cover, upward vertical 
motions on top of friction layer and net moisture convergence all re- 
vealed rather similar patterns. Figure 11 (a,b,c,d) shows the friction 
layer vertical motions in 6-hour steps with superimposed 1000 millibar 
wind directions. Results at 18 and 2b hours again appeared very realistic. 
Largest vertical motions are of the order - 300 x 10 ^ mb/sec, i.e. 2 s U cms 
rising motion. 

The evolution of the vertical motion field in a vertical plane is 
shown in figure 12. The vertical plane is a north-south cross section 
from 28°N through 28°S at 152°E longitude. This plane intersects an act- 
ive portion of the ITCZ. The field of vertical motion undergoes consider- 
able adjustment in the first 18 hours in the equatorial latitudes. This 
field appears very realistic at 18 and 2b hours of forecast time. The 
largest magnitudes of vertical motion occur at 2b hours .and are - 200 x 
10 mb/sec (* 3 cm/sec) and + 500 x 10 mb/sec (** - 7 cm/sec) in the 
mid troposphere. The adjustment of the vertical motion field produces a 
solution that is quite asimilar to initial vertical motions. Magnitudes 
of “ 10 qq are of the order of 25 x 10 - ^ mb/sec (vertical motion is zero at 
100 millibars). 

d ) Some parameters related to parameterization of cumulus convection : 

In figure 13 we show the tine variation in 6-hour steps at 152°E 
and between 28°S and 28°N of the following quantities : 

i) Specific humidity at 900 mb (gm/kgm) 

ii) Percent cloud cover 

8 1 

iii) Net moisture convergence (unit 10 mb meter sec). 
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The field of specific humidity changes very slightly during the 
first 2k hours. Along this section specific humidity varies between 10 
and 18 gm/kgm between the equatorial maximum value and lower values near 
28°N and 28°S. The region of maximum moisture is not the most convec- 
tively active region. Percent cloud cover (maximum 1 per cent) and the 
net moisture convergence are very nearly proportional to each other. 

These fields undergo rapid variations during the first 18 hours. 

It appears from the nature of these variations that spurious ver- 
tical advections of heat and moisture take place during the adjustment 
process. For short-range predictions in low latitudes, it is highly 
desirable that such spurious effects be removed right from the start. 
Unfortunately , the present model is incapable of alleviating this prob- 
lem. To the best of the author's knowledge, suppression of this initial 
oscillation has not been satisfactorily carried out in any real initial 
data experiments. Miyakoda and Moyer (1968) and Nitta and Hovermale 
(1967) have proposed procedures for initialization that could alleviate 
this difficulty, but because of limited tests in their studies, no 
general statement can be made at present. 

e) An experiment with constant relative humidity on each pressure level 
at initial time : 

As stated earlier, the initial moisture distribution is somewhat 
unreliable, hence relative humidity obtained from subjective analysis was 
averaged horizontally over the domain of integration in this experiment. 
Thus initially no gradients of relative humidity were present on constant 
pressure surfaces but it should be noted that specific humidity is not a 
constant. The results of numerical prediction of the mass and motion 
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fields at 2U hours were quite similar to those obtained from detailed 
analysis of moisture. This implies the following important feature of 
the parameterization of the cumulus scale in this experiment: 

V 3 q\V qV » W, i.e.W 3 V q 

is small. Mass convergence is a more pronounced feature of the flows 
and accounts for the net convergence of flux of moisture while hori- 
zontal advection of moisture contributes less. This statement must be 
interpreted with seme caution, because this is only true in a gross 
sense. The author has carried out another detailed synoptic study 
(fCrishnamurti , 1968b) of various parameters of large-scale flows and 
their relation to convective cloud cover as revealed by radar returns 
over the Caribbean. It was found that mass convergence and moisture 
convergence were not of the same sign over all regions near 900 mb. 

Since the moisture variable is one of the most unreliable fields in 
our input we might suggest that a constant distribution of relative 
humidity on pressure surfaces might be a useful compromise for very 
short-range prediction. 

8. Some comments on the dynamics of the intertropical convergence zone 
The field of vertical motions in the balance initial state is 
in good agreement with the synoptic disturbances north of 15°N and 
south of 15°S. In the vicinity of the equator initial vertical motion 
field shows a general region of rising motion north of the surface asymp- 
tote of convergence and does not appear to be very realistic. A short 
range prediction in a matter of 18 to 2 h hours yields a rather realistic 



pattern of rising motion that extends from east to west almost along 
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the asymptote of convergence in the surface forecast wind field. The 
initial discrepancy must be due to the inadequate dynamics of the bal- 
anced initial state which is only valid for flows characterized by a 
small Rossby number. The adjustment of the three-dimensional motion 
field to the temperature, moisture, mass field, and other external 
boundary parameters is interesting. A narrow axis of computed cloud 
cover at 2b hours along ^°N raises the question whether the output in- 
formation at 2h hours could be used to investigate some aspects of the 
dynamics of the intertropicaJ convergence zone. A short-range predic- 
tion is certainly not a simulation experiment, our initial data does 
already contain the proper ingredients in the motion, temperature and 
moisture variables to resolve in a matter of 2b hours a narrow region 
of cloud cover and rising motion along the asymptote of low-level con- 
vergence. At best, we could hope to merely carry out a diagnostic study 
with the detailed output information on tape at 2b hours . At this time 
there is some sort of a balance of the dependent variables for large 
Rossby number flows in the vicinity of the equator. 

Diagnostic studies are generally limited in scope, if calcula- 
tions are carried out in a small domain (as they are in our studies), 
any description of the local circulation can be somewhat biased. A per- 
tinent question that could be raised at this point would be on the 
maintenance of the kinetic energy of the wind field. 

Let us consider a box enclosing the region, l6°N through l6°S, 
1000 millibars through 100 millibars and from 10U°E through 128 °U. 
Because of the imposed cyclic continuity there are no boundary fluxes 



in the zonal direction. 
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The following integrals, respectively, represent the conversions 
of zonal kinetic energy and eddy available potential energy into eddy 
kinetic energy, see Lorenz (i 960 ), and Manabe and Smagorinsky ( 1967 ) • 




Equations (**5) and (1*6) represent contributions to energy transforma- 
tions at a fixed latitude, primes are departures from zonal averages 
and bar is a zonal average, other symbols are explained elsewhere. 

The net energy conversion over the box may, if desired, be obtained by 
merely integrating these conversions over the surface area of the box. 
Both of these quantities may be expressed in millibar meter sec ^ units. 
This unit seems more appropriate, because we are merely comparing the 

two rates of energy transfer with respect to each other. For similar 

2 

calculations Manage and Smagorinsky (1967) use Joule/(cm day) as the 

2 

unit, the conversion between these sets of units is 1 Joule/(cm day) = 
(g^O mb meter sec -1 . 

Figure lh shows the results of calculations at t = 2** hrs. In 
the vicinity of the intertropical convergence zone, the so-called baro- 
tropic conversion is a transfer from eddy to zonal kinetic energy while 
there is a stronger conversion of eddy potential to eddy kinetic energy. 
This result is not surprising if we examine the synoptic maps at differ- 
ent levels at t = 2U hours. 

The wind field at 1000 and 200 millibars exhibits a belt of 
easterlies near the ITCZ. Easterlies are strong at low levels and at 
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higher levels, being ouite variable and weak at the 600 and 1*00 millibar 
surfaces. At the lower level intense rising centers u' <0 are found 
where u' > 0 (weaker easterlies near the ITCZ, see figure 15) in this 
region ^ <0. At higher levels w' < 0, u' <0 and > 0. The net re- 
suit of these configurations of the motion field yields 



1 

g 



1-1000 

to'u' 

■'100 



3u 

ap 



dp < 0. 



Furthermore , along the asymptote of convergence lies a region of weak 

surface winds, figure 15. Just north of the ITCZ, along a latitude 

circle v' < 0, u' >0, ^ < 0, south of the ITCZ v' > 0, u' > 0, > 0, 

a,v ay 

and hence 



rlOOO 
g •'100 



u’ v’ 



3u 

ap 



dp < 0. 



Thus , the geometric configuration of the computed velocity field at 2h 

hours has a net barotropic energy exchange (K^ • K^,) < 0. This exchange 

is more than compensated by the ascending warm air (P^ ° K^) near the 

ITCZ, figure 11. The barotropic exchange away from the ITCZ shows 

(K^ • Kg) > 0 in both the Northern and Southern Hemispheres, this result 

is in accordance with Charney's (1963) scale analysis for large-scale 

motions in the tropics. The dominant energy conversion term is (P^ • K^,), 

though between 8°S and l6°S and between 8°N and l6°N (K • K ) is an 

L h 

appreciable fraction of the total exchange. 

The intense rising motion (u)’) along the ITCZ in our model is 
mainly attributable to the coupling of strong frictional convergence and 
the convective release of latent heat. An experiment was also performed 
without latent heat but including surface friction. The results indicated 
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a breakdown of the elongated asymptote of convergence into several large 
eddies in a matter of 18 hours. Perhaps this result was a consequence 
of improper energetics and the adjustment of our balance initial state 
to this situation. 

Because of the limited nature of the present investigation, a 
complete energetics of an open system containing the ITCZ was not carried 
out. 



10. Some concluding remarks: 

The procedure for parameterization of cumulus convection used in 
this study is rather crude. Vertical momentum transfer in sub-grid scale 
motion is not included. Recently Ooyama ( 19 67 ) and Estoque ( 1967 ) have 
considered this effect in their studies. Its inclusion would be a simple 
extension of the present study. We have shown that a somewhat realistic 
short-range prediction of larger scale flows in equatorial latitudes is 
not possible without a parameter is at ion of the cumulus scale. We have, 
however, made a tacit assumption that cumulus convection is always express- 
ible as a function of larger-scale moisture convergence. This is generally 
the case, but it excludes the possibility of isolated cumulus convection 
with local smaller-scale mass compensation. Convection over heated islands 
and sea breeze effects generally occur without bearing any relationship to 
organized larger-scale mass or moisture convergence, though they may be 
modified by the prevailing synoptic situation. This is not, however, a 
hopeless problem. Parameterization of cumulus convection could, in prin- 
ciple, be carried out to include both of these effects separately. There 
is thus an obvious need for a detailed climatological investigation of 
cloudiness over all islands and coastal regions that are known to exhibit 
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local diurnal convection. If such statistics were accumulated over the 
entire tropics then this information could be introduced in some empiri- 
cal manner into larger scale numerical prediction models. This is neces- 
sary because convective clouds, no matter how they are formed, do transfer 
large amounts of moisture, heat and momentum in the vertical direction and 
the larger-scale motions are gradually altered as a result. If a coarse 
mesh of the order of few hundred kilometers is used in prediction experi- 
ments then it may be important to include this second class of convective 
element. Parameterization of cumulus convection requires a knowledge of 
the distribution of moist unstable regions in the horizontal and in the 
vertical. A knowledge of the distribution of the moisture variable is 
thus quite important, yet we find that the initial analysis of this field 
is very difficult especially from the currently available data networks. 
Several massive research efforts in the tropics are currently being con- 
templated in the United States and other countries. It would be highly 
desirable to examine ways to obtain detailed three-dimensional distribu- 
tion of moisture variables in low latitudes. 

Some of the important results of the present study are as follows: 

The dynamical structure of the tropical atmosphere in equatorial 
latitudes is poorly described by the balance model, but is reasonably well 
represented if a short-term numerical prediction is carried out. 

At 2h hours of forecast time the following features are of inter- 
est. A well marked surface asymptote of convergence, the ITCZ, is found 
at h°N. Along this asymptote an axis of frictional mass and moisture con- 
vergence is found with a maximum computed convective cloud cover. This 
study reveals a strong coupling between these three features. Along the 
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equator, sinking motion and absence of clouds is indicated by the compu- 
tations. Furthermore we find that in the vicinity of the ITCZ, energy 
transformations show that (P^, • K^,) is positive, in agreement with the 
results of Fanabe and Smagorinsky (1967) but we also find a surprising 
result that (Kg • K^,) is positive. We have shown that the latter result 
is possible because of the detailed structure of wind field in low lati- 
tudes. Away from the ITCZ we find that (K„ » K ) is positive implying a 
barotropic instability of the zonal current. This study has demonstrated 
a limited success in short-range numerical prediction from real initial 
data in low latitudes. 

Besides the difficulties of the parameterization of the cumulus 
scale other major defects of this study are related to lack of observa- 
tions, contamination of forecasts due to lateral boundary conditions and 
the presence of gravity inertio-oscillations during the ^initial phase of 
the forecast period. Further work is continuing to examine in detail 
some of these problems. 



A c kno wl e dg eme nt s 

Mr. Edward B. Rodgers, a graduate student at Florida State University, 
provided considerable assistance in data processing and analysis. Initial 
financial support for this work was provided by Air Force Cambridge Research 
Laboratories, Grant No. MIPR ES-7-967, while the author was at Naval Post- 
graduate School, Monterey, California. This work is being continued under 
support by the Atmospheric Science Section, National Science Foundation, 

NSF Grant GA-1U80, at Florida State University. Computations were carried 
out on the IBM 709^ computer at the University of California and on the CDC 
6U00 computer at Florida State University. 



REFERENCES 



Baumhefner, D. P. 

Benwell , G . R . R . and 
Bretherton, F. P. 

Berkofsky, L. and 
Bertoni, E. A. 

Bushby, F. H. and 
Timpson, M. S. 

Charney, J. 

Charney, J. G. and 
Eliassen, A. 

Dean, G. A. 

Elsasser, H. W. 
Estoque, M. A. 

Grammeltvedt , A. 

Gray, W. M. 



1968 'Application of a diagnostic Numerical 
Model to the Tropical Atmosphere , ' 

Mon . Wea . Rev . , £6, pp. 218-228. 

1968 'A pressure oscillation in a 10-level 

atmospheric model , ' Quart. J. R. Met . 
Soc . , 9k, pp. 123-131. 



i 960 'Topographic charts at one degree inter- 
section for the entire earth,' G.R.D . 
Research Notes. No. .42 , pp. 1-^3, Air 
Force Cambridge Research Laboratories, 
Bedford, Mass. 

1967 'A ten-level atmospheric model and frontal 

rain,' Quart. J. R. Met. Soc ., 93 , 
pp. 1-17 • 

1955 'The use of primitive equations of motion 

in numerical prediction,' Tellus , 7_, 

pp. 22 - 26 . 

1963 'A note on large-scale motions in the tropics,' 
J. Atmos. Sci . , 20, pp. 607-609. 

I 96 U 'On the growth of the hurricane depression,' 

J . Atmos . Sci . , 21 , pp. 68-75- 

1956 'The 1955 mean monthly wind circulation over 

the tropical central Pacific area,' 

Institute of Geophysics, University of 
California , lb pp. 

1968 ’Diagnosis of early baroclinic NWP models,' 

J. Appl. Met . , 7_, pp. 153-159- 

1967 'Exchange processes due to penetrative 

convection,' Bull. Araer. Met. Soc . , U8, 
p. 62 U. 

1968 'A survey of finite difference schemes for 

the primitive equations for a barotropic 
fluid , ' Technical Report, National Center 
for Atmospheric Research , Boulder, Colorado. 

1967 'Global view of the origin of tropical dis- 
turbances and storms , ' Atmospheric Science 
Paper No. llU , Colorado State University , 
Fort Collins, Colorado, 



1 * 1 * 



Hawkins , H . E . and 
Rosenthal, S. L. 


1965 


'On the computation of the streamfunction 
from the wind field,' Mon. Wea. Rev., 
93, pp. 21*5-252. 


Hinklemann, K. 


1951 


'Der Mechanismus des meteorologischen Larmes,' 
Tellus, 3, pp. 285-296. 


Jacobs, W. C. 


1951 


'Large scale aspects of energy transformations 
over the oceans,' Compendium of Meteorology 
American Meteorological Society, Boston, 
Mass. 


Krishnamurti , T. N. 


1962 


Numerical integration of primitive equations 
by a quasi-lagrangian advective scheme, 1 
J. Appl. Met., 1, pp. 508-521. 




1968a 


’A diagnostic balance model for studies of 
weather systems of low and high latitudes , 
Rossby number less than l, f Mon* Wea. Rev., 

96 , P p. 197 - 207 . 




1968b 


'A calculation of the percent area covered 
by convective clouds from moisture con- 
vergence,' J. Appl. Met., 7, PP. 181*-195» 


Kuo, H. L. 


1961 


'Convection in conditionally unstable atmos- 
phere,' Tellus, 13, pp. 1*1*1-1*59. 




1965 


'On formation and intensification of tropical 
cyclones through latent heat release by 
cumulus convection,' J. Atm. Sci. , 22, 
pp. 1*0-63. 


Leith, C. 


196U 


’Lagrangian advection in an atmospheric model,' 
Tech. Note No. 66, pp. 168-17 6, World Met. 
Organiz., Geneva, Switzerland. 


Lilly, D. K. 


1965 


'On the computational stability of numerical 
solutions of time dependent non-linear 
geophysical fluid dynamics problems,' Mon. 
Wea. Rev., 93, pp. 11-26. 


Lorenz, E. N. 


I960 


’Energy and numerical weather prediction,’ 
Tellus, 12, pp. 361*— 373. 


Manabe , S . and 
Smagorinsky, J. 


1967 


'Simulated climatology of a general circula- 
tion model with a hydrologic cycle II: 
Analysis of the tropical atmosphere,' 

Mon. Wea. Rev. , 95, pp. 155-169* 


Matsuno, T. 


1966 a 


'Quasi geostrophic motion in the equatorial 
area,' J. Met. Soc. Japan., i+U , pp. 25-1*3- 



Matsuno, T. 



Mintz, Y. 



Miyakoda, K. and 
Moyer, R. W. 



Miyakoda, K. , 

Smagor insky J. , 
Strickles, R. F 
Hembree, G. D. 

Morskoi Atlas 

Nitta T. and 

Hoverwale, J. B 



0kland, H. 



Ooyama, K. 



Phillips, N. A. 



Riehl, H. 



1966b 

196k 

1968 
1968 

. and 

1953 
1967 

1965 
1963 

1967 

1956 
I960 

191*5 



'A finite difference scheme for time integra- 
tions of oscillatory equations with second 
order accuracy and sharp cut-off for high 
frequencies , ' J . Met . Soc . J apan , M*_, 

pp. 85-88. 

'Very long-term global integration of the 

primitive equations of atmospheric motion,' 
Techn i cal Note No. 66 , pp. lUl-l67, World 
Meteorological Organization , Geneva, 
Switzerland. 

'A method of initialization for dynamical 
weather forecasting,' Tellus , 20 , pp. 113- 
128. 

'Experimental extended predictions with a 
nine level hemispheric model. ' Technical 
Report, U.S. Department of Commerce , pp. 
1-53, ESSA, W a shington, D. C . 

Russian Marine Atlas , Vol. II , Plates I-85. 

'On analysis and initialization for the 
primitive forecast equations,' Weather 
Bureau Technical Memorandum No. EMC, h2 , 
pp. 1-38. National Meteorological Center, 
Suitland, Maryland. 

'An experiment in cyclogenesis prediction by 
a two-level model , Mon. Wea. Rev . , 93 , 

pp. 663-672. 

'A dynamical model for the study of tropical 
cyclone development , ' Dept, of Meteorology 
and Oceanography, New York University , 

26 pp. 

'Numerical simulation of the life cycle of 
tropical cyclones , ' NSF Report No. GA-623 , 
New York University , University Height, 

New York. 

'The general circulation of the atmosphere, 
a numerical experiment , ' Quart . J . R . 

Met . Soc . , 82, pp. 123-161* . 

'On the problem of initial data for primitive 
equations,' Tellus , 12, pp. 121-126. 

Tropical Meteorology , McGraw-Hill Book Co., 

New York. 



Rohwer, C. 



Rosenthal, S. L. 



Sadler , J . C . 



Shuman, F. G. and 
Hovermale , J . B . 

Smagorinsky , J. and 
Staff Members 



U.S. Bureau of 
Commercial 
Fisheries 



Webb, K. E. 



1931 'Evaporation from free water surfaces,' 

U.S. Dept. Agriculture, Technical Bulletin 
No. 217 , 96 pp. , Washington. 

1965 'Some preliminary theoretical considerations 

of tropospheric wave motions in equatorial 
latitudes,' Mon. We a. Rev . , 93 , pp. 605- 
612. 

1967 'The tropical upper tropospheric trough as 

a secondary source of typhoons and a 
primary source of tradewind disturbances,' 
Final Report , No. AFCRL-67-0203 , pp. 1-U8, 
University of Hawaii, Institute of Geophys . 
to AFCRL, Bedford, Mass. 

1968 'An operational six layer primitive equation 

model,' J. Appl. Met . , 7_, (to be published) 

1965 'Prediction experiments with a general circu- 
lation model , ' Proceedings of the Interna- 
tional Symposium on Dynamics od Large Scale 
Processes in the Atmosphere , pp. 23-30, 
Moscow, USSR. 

1962 'Monthly mean charts, sea surface temperature. 
North Pacific,' U.S., Fish Wildlife Service , 
Circular No. 13^, Biological Laboratory , 
Stanford, California. 

1960 'An inverstigation of the evaporation from 

Lake Eucumbene , ' Technical Paper No. 10 . 
1-75 pp. Commonwealth Scientific and 
Industrial Research Organization, Australia 
Division of Meteorology and Physics, 
Melbourne, Australia. 

1961 'A detailed analysis of typhoon formation, 

J. Met. Soc. Japan , 39 , pp. 187-21 k. 



Yanai , M. 



I 



Figure Captions 



Figure 1. Sketch of nine points on a pressure surface used in quasi- 
Langrangian advection. P is the origin of a parcel at time 
t - At, and Q is a grid point where parcel arrives at time t. 

Figure 2. Vertical staggeriing of variables in the initialization and 
prediction phase. 



Figure 3. Procedure used for the parameterization of cumulus convection. 



Figure U. (a) (top) Streamlines (solid lines) and isotachs (dashed 

lines , knots ) for March 1 , 1965 , 1200Z , at 
1000 millibars, t = 0 in our experiment. 

(b) (bottom) Streamlines (solid lines) and isotachs (dashed 
lines, knots) for March 1, 1965, 1200Z, at 200 
millibars, t = 0 in our experiment. 



Figure 5* (a) (top) Specific humidity distribution at 900 millibars, 

units gm/kgm, for March 1, 1965, 1200Z, t = 0 in 
our experiment. 

(b) (bottom) Temperature distribution at 500 millibars, 

units °A, for March 1, 1965, 1200Z, t = 0 in our 
experiment. 



Figure 6. (a) (top) Satellite mosaics for March 1, 2 and 3, 1965, in 

the vicinity of lUo°E longitude. 

(b) (bottom) Satellite mosaics for March 1, 2 and 3, 1965, in 
the vicinity of 1T0°E longitude. 
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Figure 7* (a) (top) Initial streamfunction at 200 millibar surface, 

units 10^ meter^ sec ^ , March 1, 1965, 1200Z. 

(b) (bottom) Initial stream function at 1000-millibar surface, 
units 10^ meter^ sec March 1, 1965, 1200Z. 

Figure 8. Initial balance vertical motions (solid lines) at 900 millibar 
units 10^ millibar/sec. The initial calculated percent convec- 
tive cloud cover is shown by dashed lines and this region is 
shaded. March 1, 1965, 1200Z. 

Figure 9- (a) (top) 24-hour forecast of the wind field at 1000 milli- 

bars, solid lines are streamlines and dashed 
lines are isotachs in knots. The plotted winds 
are observations at verification time, i.e. March 
2, 1965, 1200Z. 

(b) (bottom) Same as figure 9(a) except at 200 millibars. 

Figure 10. (a, b, c and d) 6-, 12- , 18- and 24-hour forecast of the dis- 
tribution of percent area of convective cloud 
cover, shaded region. The wind arrows are the 
directions of the forecast wind vectors at the 
1000 mb surface. 

Figure 11. (a, b, c and d) 6-, 12- , 18 - and 24-hour forecast of the dis- 
tributions of 900-millibar vertical motions, 
units 10~5 millibar/sec, the shaded region indi- 
cates rising motions. The wind arrows are the 
directions of the forecast wind vectors at the 



1000 mb surface. 



Fig. 12. 



Figure 13 



Figure l 1 * 



Figure 15 



Vertical cross-section of vertical motion (10 ^ millibar/sec) 
at 6, 12, 18 and 2k hours of forecast time. This cross 
section is taken in a south-north plane at 152°W longitude, 
and it intersects the ITCZ. 

Time variation (along ordinate) of selected parameters at 
152°W from south to north (along abscissa). 

q specific humidity at 900 mb gm/kgm 
a percent cloud cover 

8 l 

I net moisture convergence (10 mb meter sec). 

Selected energy conversions at t = 2k hours of forecast time, 
evaluated from equations (1*5) arid (1*6) of text. 

An enlarged view of the 1000 mb, 2U-hour forecast of wind direc- 
tions( arrows) and wind speeds (knots, solid lines) in the 
vicinity of the ITCZ (heavy dark line). 
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Figure 1. Sketch of nine points on a pressure surface used in quasi- 
Langrangian advection. P is the origin of a parcel at time 
t - At , and Q is a grid point where parcel arrives at time t. 
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Figure 2. Vertical staggeriing of variables in the initialization and 



prediction phase. 



FLOW CHART FOR CUMULUS SCALE HEATING:- 




Figure 3* Procedure us^d for the parameterization of cumulus convection. 




Figure 4. (a) (top) Streamlines (solid lines) and isotachs (dashed 

lines, knots) for March 1, 1965, 1200Z, at 
1000 millibars, t = 0 in our experiment. 




(b) (bottom) Streamlines (solid lines) and isotachs (dashed 
linos, knots) for March 1, 1965, 1200Z, at 200 



millibars , t = 0 in our experiment 





our experiment. 




units °A, for March 1, 1965, 1200Z, t = 0 in our 



experiment. 
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Figure 6. (b) Satellite mosaics for March 1, 2, and 3, 1965 in vicinity of 170°E longitude. 




Figure 7. (a) (top) 



Initial streamfunction at 200 millibar surface, 
units 10^ meter^ sec \ March 1, 1965, 1200Z. 




(b) (bottom) Initial stream function at 1000-millibar surface, 

5 2-1 

units 10 meter sec , 



March 1, 1965, 1200Z 




shaded. March 1, 1965, 1200Z. 
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(top) 24-hour forecast of the wind field at 1000 millibars, 
solid lines are streamliner and dashed lines are 
isotachs in knots* The plotted winds are observations 
at verification time, i.e. March 2, 196 5, 1200Z. 
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(b) (bottom) Same as figure 9(a) except at 200 millibars. 
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Figure 10« l •• , v >> n * d a) 6-, 12, 18, and 24-hour forecast of the dis- 



cribut io.i of oercent area of convective cloud cover, 
shaded region* The wind arrows are the directions of the 




t wind vectors at the 1000 mb surface 
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Figure 10* (d) 
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Figure li. 



(a, b, c, & d) 6-, i2 T , 18-^ ana 24-hour forecast ofj.the dis- 
tribution of 900 -millibar vertical motions , units lO’“' > millibar/ 
sec, the shaded region indicates rising motions* The wind 
arrows are the directions of the forecast wind vectors at the 
1000 mb surface. 
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F igure 1 1 * ( c ) 
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Fig. 12. Vertical cross-section of vertical motion {10 - ^ millibar/sec) 
at 6, 12, 18 and 2U hours of forecast time. This cross 
section is taken in a south-north plane at 152°W longitude, 
and it intersects the ITCZ. 
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Figure 13* Tim* variation (along ordinate) of selected parameter* »t 
152 V froa south to north (along abscissa). 
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Figure lU. Selected energy conversions at t = 2k hours of forecast time, 
evaluated from equations (45) add (46) of text. 
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